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ABSTRACT 

We discuss high-resolution simulations of one-dimensional, plane-parallel shock 
waves with mean speeds between 150 and 240 kms~^ propagating into gas with Alfven 
velocities up to 40 kms~^ and outline the conditions under which these radiative shocks 
experience an oscillatory instability in the cooling length, shock velocity, and position 
of the shock front. We investigate two forms of postshock cooling: a truncated single 
power law and a more realistic piecewise power law. The degree of nonlinearity of the 
instability depends strongly on the cooling power law and the Alfven Mach number: for 
power-law indices q < typical magnetic field strengths may be insufficient either to 
stabilize the fundamental oscillatory mode or to prevent the oscillations from reaching 
nonlinear amplitudes. 



Subject headings: hydrodynamics — instabilities — shock waves 



Introduction 



Work by Langer, Chanmugam, & Shaviv (1981), Chevalier & Imamura (1982), and Imamura, 



Wolff, & Durisen (1984) shows that the physical structure of some radiative shocks and their 
cooling regions is subject to a cooling overstability. The shock velocity and the length of the 
cooling column of postshock gas may oscillate instead of evolving in a steady manner. The linear 
perturbative analysis of Chevalier &; Imamura, which assumes power-law cooling laws A cx p^T", 
indicates that one-dimensional nonmagnetic shocks are stable when a is sufficiently large, a ^ 0.8. 
If the cooling rises rapidly with temperature, perturbations to a cooling region in steady state 
that increase the shock velocity Vs cause more rapid cooling and require shorter cooling lengths, 
and the shock is stable against perturbations. On the other hand, if the cooling law does not rise 
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sufficiently rapidly with temperature, then the cooling length oscillates overstably between fast 
shocks with long cooling lengths and slow shocks with short cooling lengths. 



Toth & Draine ( 1993| ) extend the perturbative analysis to include a frozen one-dimensional 



magnetic field, oriented perpendicular to the gas motion and uniform in the upstream gas. 
Compression at the shock front tends to align the postshock magnetic field with the plane of the 
shock. The field's contribution to the total pressure tends to stabilize oscillations. 

In this paper we discuss sets of one-dimensional, plane-parallel simulations of radiative 
shocks, both with and without the transverse magnetic field. Shocks in the interstellar medium 
propagating faster than approximately 150kms~^ are candidates for these instabilities. For 
example, the blast waves of supernova remnants may pass through a velocity range in which 
the shock velocities are high enough and the cooling times short enough, compared with the 
remnant age, for these oscillations to occur ( [Kimoto fc Chernoff 1997| ). For cooling laws we adopt 
A = n^/(T), where rin is the hydrogen- nuclei density. For the cooling function f{T) we use either 
(1) single power laws or (2) a piecewise power-law fit (for temperatures above 3 x 10*^ K) to the 
result of Raymond, Cox, & Smith ( |1976D . 

We note that the latter cooling function assumes that the emitting plasma is in collisional 
equilibrium. However, because collisional, ionization, and recombination rates are slow compared 
to the cooling rate, an accurate evaluation of the cooling necessitates a time-dependent treatment 
of the ionization evolution. Several specific cases of one-dimensional simulations incorporating 
explicit evolution of ionization states are discussed by Innes, Giddings, & Falle ( |1987| ), Gaetz, 
Edgar, & Chevalier (|l988D , and Innes ( ^99^ ). Of these three papers, only Innes (|1992|) incorporates 



a transverse magnetic field. Such treatments are computationally very demanding. By using 
simpler cooling laws, we are able to investigate a wider range of parameter space in mean shock 
speed and magnetic field strength, as well as assess the sensitivity of the results to the form of the 
cooling law. Our qualitative conclusions should apply to more accurate treatments as well. 



Our piecewise power-law cooling function differs from that of Toth & Draine ( |1993 ) as 



discussed below, and our results differ qualitatively from their similar suite of calculations. 
Because our cooling function is more unstable for mean shock velocities between 150kms~^ and 
210kms~^, we find that the stabilizing effect of magnetic fields is reduced. In our analysis, larger 
magnetic fields are necessary to stabilize large-amplitude, fundamental-mode oscillations, and as 
well to stabilize all oscillatory modes. We do not claim that our treatment leads to results that 
are physically more reliable: the differences between the two sets of results serve to point out that 
quantitative results will require a more precise treatment of the cooling. 

Most recently, Walder Sz Folini ( |1996| ) discuss spherically symmetric, one-dimensional, 
high-resolution simulations of the radiative cooling instability. Like ours, these simulations use a 
piecewise power-law cooling function. They identify five qualitatively distinct oscillation types. 
We note that the two types that we discuss in some detail below correspond to the "strong forms" 
that they describe for the fundamental and first overtone modes. We do not draw further parallels 
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with their work, however, since their simulations do not include magnetic fields, an essential 
ingredient in our discussion. 

In the following section we discuss simulations carried out with our two cooling functions 
in parameter spaces in which we vary the magnetic field strength and properties of the cooling 
function. We identify two distinct types of nonlinear oscillations: one in which the fundamental 
mode dominates the motion, and one in which the first-overtone mode dominates. In the 
final section we summarize and draw conclusions based on our results, taking into account the 
limitations of the systems we have studied. 



2. Simulations of oscillations 

2.1. Equations of motion and initial conditions 

In order to study the oscillatory cycles of this cooling over stability, we consider one- 
dimensional, plane-parallel systems of ideal (inviscid and with zero thermal conductivity), 
infinitely conducting gases, with a magnetic field perpendicular to the direction of motion. Our 
cooling law takes the form A = n^/(T). To simplify the investigation of the oscillatory behavior, 
in some cases we adopt for f{T) single power laws (with low-temperature cutoffs). In the other 
cases, as an approximation to a realistic cooling function, we use for f{T) a piecewise power-law 
fit to the result of Raymond et al. ( 1976| ). To account for the turnoff of cooling when the gas 



reaches a sufficiently low temperature to recombine, we arrange for the cooling to vanish below a 
cutoff temperature Tc ~ 2 x 10'' K. Figure || shows the function we adopt, as well as the function 
used for application to the interstellar medium by Toth & Draine ( 1993| ). The points indicated in 



the figure denote the shock temperatures of steady-state shocks (in the absence of magnetic fields) 
with velocities 150, 180, 210, and 240 kms"^ 

Since the cooling law takes this simple functional form, the length and time scales of solutions 
are simply proportional to 1/?t-h, and all densities, pressures, and energies proportional to nj|. 
Because of this scaling, the particular value we choose for the upstream density is arbitrary, 
and so usually we quote values in terms of nJJ (omitting the superscript "in" where there is no 
ambiguity) . 

Our numerical calculations use an Eulerian finite-difference algorithm that can evolve flows 
with strong shocks and gradients. We use operator splitting to separate the cooling from the 
magnetohydro dynamics, in which we evolve four conserved densities {p, pv, etoti and B) with 
a flux-corrected transport scheme ( ^alesak 1979| ). A multigrid method using some ideas of 
Berger & Colella ( |1989| ) allows us to place high resolution on localized parts of the flow as required 
for numerical accuracy. Details of the methods may be found in the Appendix of Kimoto & 
Chernoff ( |l99^ ). 

In all of the simulations, supersonic gas of temperature 2500 K flows uniformly into one end 
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of the computational domain and is shocked. This gas is assumed to be completely preionized by 
the photons emitted near the shock and to contain helium in a 1:10 ratio to hydrogen. (In our 
simulations, we give the inflowing gas a hydrogen- nuclei density of = 50 cm~^, but because 
of the scaling properties of our equations of motion, the particular value has no significance.) 
Although the temperature of our unshocked gas is unrealistic both for preionized gas and for a 
quiescent interstellar medium, the dynamics of the shock and the cooling column is unaffected as 
long as the shock is strong. The gas cools and accumulates near a wall at the other end, where 
we impose a perfectly reflecting boundary condition. When the cooling region evolves in a steady 
state, the quiescent, fully cooled gas builds up at the wall, and the cooling region moves outward 
at constant velocity. We find considerably more complicated dynamics. 

The initial conditions are constructed to correspond closely with the steady-state cooling 
flow for a given shock speed and upstream Alfven velocity. The main ingredient is a numerical 
solution to the equations of motion, generated by working in the shock frame and assuming no 
time dependence. First, this steady-state solution is truncated at a point where the gas has 
approximately attained its final cooled, dense state (beyond the main part of the cooling column), 
and second, the transition between upstream and shocked gas is smoothed slightly. The numerical 
evolution scheme smooths discontinuities over several gridpoints. By necessity the smoothing that 
we employ in generating initial conditions differs slightly from the smoothing developed under 
the numerical scheme. We then transform the solution into the frame of the cold, dense gas that 
builds up at the wall boundary. In most cases the inflow velocity is only slightly less than the 
average shock speed since the gas compresses strongly. During the highly nonlinear oscillations 
the conditions of the cooling column do not resemble those of the steady state, and waves are 
driven downstream into the cooled, dense gas. We find, however, that the measured average 
shock propagation speed fluctuates only slightly about the shock speed of the steady-state initial 
condition. Table |l| lists the results for twenty cases discussed below — covering mean shock speeds 
from 150 to 240kms~^ and upstream Alfven velocities up to 20kms~^ — that use the piecewise 
power-law cooling function. 



2.2. Parameter space and instability types 

The linear perturbative analysis of Chevalier & Imamura ( |1982| ), which assumes f{T) oc T" 
and no magnetic field, yields a discrete spectrum of modes for motions of the cooling column. The 
fundamental mode has frequency 

I'F ~ 17— , (1) 

and the overtone frequencies are approximately odd multiples of the fundamental. The overtone 
modes are unstable over a greater range of cooling power laws than is the fundamental. 

A transverse magnetic field can stabilize the oscillations. The perturbative analysis of Toth & 
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Draine ( |1993| ), which incorporates such a magnetic field, uses a coohng function 

/td(T) OC (2) 



for some turnover temperature Tf. For shock temperatures much greater than Tt, the minimum 
stabilizing magnetic field depends mainly on a: the field is greater for smaller values of a (i.e., 
the more unstable cases), and it is greater for overtone modes than for the fundamental. In 
addition, the stabilizing field for all overtone modes is approximately the same. In nonlinear 
calculations typically only the first overtone mode is observed. The results can be summarized 
as follows: the minimum stabilizing field for mode i may be parameterized by Mj^^ (where 
the Alfven Mach number is Ma = Vs/va, the Alfven velocity is va = B /(ATTpi^^)^/'^ , and the 
upstream density is pi^). For small magnetic fields, MJ^ < Mj^p < M^]j, both fundamental and 
overtone modes are unstable; for intermediate fields, MJ^p < MJ^ < MJ^q, only the overtones 
are unstable; and for large fields, M^^p < M^]j < M^^ , all modes are stable. The short dotted 
lines in Figure |2| show approximately the Mj^ (lower) and MJ^q (upper) curves in the (a, M^^) 
parameter space as calculated by Toth & Draine using the cooling function (eq. [Q) with turnover 
temperature Tt = W~'^fiVg/kB, where fi is the mean particle mass and Vs is the mean shock speed. 



2.3. Oscillations without magnetic fields 

As background, we begin by describing properties of oscillations in the absence of magnetic 
fields. For simulations adopting the piecewise power-law cooling function. Table |2| lists typical 
amplitudes (minimum to maximum shock velocities) and periods for a range of shock speeds — from 
130 to 240 km s"^— that should cover a range of supernova remnants with oscillating shocks 



(Kimoto & Chernoff 1997). The corresponding range of shock temperatures covers a break in 
the power-law fit (locally proportional to T") that enters our cooling function: Above 5 x 10^ K, 
which corresponds to a 190 kms~^ shock (neglecting magnetic fields), a ~ —0.1; from 2.5 x 10^ to 
5 X 10^ K, we adopt a ~ —2.2 (for which the oscillations are more unstable). 

From these simulations we make several observations that are generally consistent with 



the results of previous work done with a variety of steep cooling laws (e.g., Imamura, Wolff, & 


Durisen 1984|; 


Innes, Giddings, & Falle 1987; 


Gaetz, Edgar & Chevaher 1988 


; Walder & 



Folini 199^ ). The instability is so strong that within a few periods the oscillations develop and 



their amplitudes quickly saturate. When the oscillations become nonlinear, the fundamental 
mode dominates. Our amplitudes and periods typically vary by up to 10%, but the values 
indicated in Table |2| should be representative. The period of the fundamental mode is close to 
Tp = 21Lcooi/ {vs), approximately the value (cf. eq. |j|]) predicted by linear analysis ( [Chevalier & 



Imamura 1982D . We show this approximate value for comparison in Table |2|, where (in the spirit 



of small-amplitude analysis) we take Lcooi from the steady-state solutions. (In the perturbative 
treatments the gas cools to T = [for analytic convenience]; our cooling function turns off in 
a tapered manner at a finite temperature [in part, for numerical convenience]. Because the gas 
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cools rapidly as it reaches high densities, it is still possible to identify a "cooling length" without 
much ambiguity.) In our simulations the oscillations become highly nonlinear, and the oscillatory 
periods are somewhat longer than the linear-analysis values. 



2.4. Oscillations affected by magnetic fields 

Next we consider simulations with a transverse magnetic field. For simplicity we begin by 
considering cooling functions with single power laws in order to investigate the effects of power-law 
slope and magnetic field on the selection of the dominant mode. For a cooling function (cf. eq. Q) 
we use 

rr" \iT>Tt, 

^PL^^^^^jTV^Ii + itanhKr-^TM^r,)]} ifT<r,. (3) 
The quantity in braces arranges for the cooling to effectively vanish at small temperatures 



(r < 9Tt/14). For a > —0.5 we (like Toth &: Draine | 1993 |) set the turnover temperature so 
that kBTt/{iivl) = 10^3 Yot a < -0.8, however, we use ksTt/ {iiv^g) = 10 ^ for numerical 
convenience. 

In the previous subsection we noted that the fundamental mode dominates the oscillations in 
the absence of magnetic fields {M^^ =0). At several values of the power-law index a (namely, 
— 1, —0.8, —0.5, —0.3, 0, and 0.3), we bracket the reciprocal Alfven Mach number at which 

dominance switches from the fundamental to an overtone mode. Simulations yield the limits in 

shown as line segments in Figure |2|, and the curve in long-dashed lines (drawn only for 
illustrative purposes) shows approximately the behavior of M.~a ■ Naturally this curve must lie 
below the lower short-dashed line, which denotes MJ^^, the value at which the fundamental mode 
is marginally stable in the linear analysis of Toth & Draine. We note that the magnitude of ■M.'/^ 
curve is clearly distinct from that of M^^p. 

Toth &; Draine discuss the results of hydrodynamical simulations with the values of a > 
and M^-^ > 0.03 denoted by dots in Fi gure 0. They observe that the first-overtone mode 
dominates whenever the flow is unstable and that the oscillations have amplitudes ^ 10%. We 
have performed several simulations in the vicinity of the dots using our version of the power-law 
cooling function (eq. [^), and these confirm their conclusions. However, Figure ^ also shows 
that the fundamental mode dominates the oscillations over an astrophysically significant range of 
magnetic fields when a < 0. Even the simplest forms of cooling functions (e.g., as shown in Fig. |l|) 
require temperature ranges in which locally a < 0. In addition, for a < the oscillations are 
decidedly nonlinear: when the fundamental mode dominates, the velocity fluctuations (minimum 
to maximum) are at least 50% of the mean shock velocity. In the handful of cases we have 
simulated in which the overtone mode dominates (i.e., with MJ^ slightly greater than A4^^), 
the velocity fluctuations are typically ~25%. Figure ^ shows the magnitude of fluctuations in 
cooling-column length and shock velocity (relative to the steady-state values) as a function of 
when a = —0.8. When the fundamental mode dominates, as occurs for a ^ 0.08, the fluctuations 
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are very strong. The amplitudes drop abruptly when the overtone mode becomes the dommant 
mode (at M^^ « 0.1) All our numerical evidence confirms that the flows have strong, nonlinear 
fluctuations when the physical parameters lie in the region below Ai^^ in Figure |2[ Above A4^^ 
the fluctuation amplitudes, while smaller than below, may still be appreciable. The parameter 
space covered by the simulations of Toth & Draine does not cover the full range of plausible 
cooling-law slopes, and consequently several conclusions for our more realistic cooling function 
contradict their general results. 

We end our discussion of single power-law cooling functions with comments about the 
robustness of the values of with respect to changes in simulation parameters. As noted 

above, in the cases with a < —0.8, for numerical convenience we use a value of the turnover 
temperature Tt ten times larger than the higher-a cases. We have verified in the case of a = —0.5 
that the selection of the dominant mode discussed below is insensitive to this increase in T^. 
For smaller values of a, we are interested in the behavior at increasingly large magnetic fields. 
The effect on the dynamics of the cooling column should be negligible since Tj sets the thermal 
pressure of the cooled gas, but near the mode transition this thermal pressure is dominated by the 
magnetic pressure. For larger values a ^ 0, cases with small or zero magnetic fields are of interest, 
and the choice of final gas temperature does qualitatively affect the oscillations. 

We believe that our choice of boundary conditions has little effect on the value of ■ 
However, simulations of one-dimensional cooling shocks without magnetic fields by Strickland &: 
Blondin ( |1995| ) suggest that boundary conditions may affect the stability or the dominant unstable 
mode. They offer as a possible explanation that a reflecting wall may return waves to the cooling 
region and interfere with the growth of unstable modes. We use the case of q = —0.5 to check 
the dependence of ■M.'^^ on boundary conditions. In order to remove the effect of reflected waves, 
we separate the cooling region from the wall with a column of quiescent gas large enough so that 
waves generated in the region of the shock are not reflected back to the cooling column during the 
simulation. We find no change in the qualitative behavior of cases on either side of our inferred 
M.^^ boundary. 

Next we return to the more complex piecewise power-law cooling function shown in Figure |l[ 
In Figure § we show the shock-velocity fluctuations for the cases mean shock speeds 150, 180, 210, 
and 240 kms~^, each with Alfven velocities 0, 2, 5, 10, and 20 kms^^. The graphs show the shock 
velocity (in the frame of reference of the reflecting wall) as a function of time. Over this range of 
shock speeds, all -B = flows are unstable. Since the local slope of the cooling function varies, 
the lower shock speeds should be more unstable than the higher ones. We discuss the effects on 
stability, dominant mode, and oscillation amplitude as the magnetic held is increased. 

Inspection of Figure ^ shows that (1) the fundamental mode dominates over most of the 
parameter space, (2) the amplitude of oscillations is nonlinear (~ 100%), and (3) the magnetic 
flelds qualitatively alters the instability only at the lowest and highest shock speeds. In a few 
cases starting from near steady-state conditions, the growth rate of an unstable fundamental mode 



- 8 - 



is relatively slow, and the qualitative long-term behavior sets in only after 10-20 cyclic times. 
Before this occurs, the shock velocity oscillates with ever-increasing amplitude, and the variations 
in shock velocity are dominated by the overtone modes. We now discuss these observations in 
light of the conclusions drawn from the single power-law cooling simulations described above. 

In two high-speed cases [vg = 240kms~^, va = 10 and 20 km s~^), the magnetic field is 
sufficiently strong that the fundamental, but not the overtone, mode is stabilized. Under these 
circumstances the oscillations are noticeably smaller in amplitude than in the corresponding cases 
with lower magnetic fields. As expected for the first-overtone mode, the oscillation frequency 
is roughly three times that of the lower-magnetic-field 240 km s^^ cases, whose fundamental 
modes are unstable. (At high Alfven Mach numbers the contribution of the magnetic field 
to the total pressure is small in gas that has not yet become cool and dense, and the cooling 
length and mode frequencies are not much altered from the f a = case.) The simulations show 
that 0.02 < M^^ < 0.04. 

We may compare this with the results obtained using single power-law cooling. The relevant 
parameter for comparison is the power-law index a. At the shock temperature corresponding 
to 240 kms^^, our cooling function locally has a ~ —0.1. In that case, our power-law results 
(cf. Fig. 1) suggest that M^^ « 0.02, roughly consistent with the 240 km s ^ simulations. The 
value of .M.^^ is likely increased because the cooling function is steeper (a ~ —2.2) in a wide range 
of temperatures below the shock temperature. 

In one case {vg = 150kms~^, va = 20kms~^) in Figure the instability is absent; the 
oscillations shown are extremely small in amplitude and decaying. The simulations show that 
0.07 < {M^^p, M^]j} < 0.13. At a shock speed of 150kms~^, the shock temperature — if 
B = — implies a local slope a ~ —2.2. Although such a steep power law is not considered by 
the linear analysis of Toth & Draine, extrapolation from their results indicates that M^^ > 0.4 
would be required to stabilize the fundamental mode. In this case the observed stabilization is of a 
different nature: the upstream magnetic field is sufficiently large to reduce the shock temperature 
by more than 10% as compared with the B = case. At this reduced temperature the cooling 
function is shallow {a = 0) and less unstable. The combination of the curvature of the cooling 
function and the lowered shock temperature stabilizes fundamental and overtone modes in part of 
parameter space. 

Finally, we come to the issue of stability over the large central range of shock speeds. We 
see no evidence of stabilization at Vg = 180 and 210kms~^ for va < 20kms~^ in Figure |. We 
have performed additional simulations with larger magnetic fields {va = 30 and 40kms^^), and 
the results for the fundamental mode are summarized in Figure The vertical line segments 
indicate constraints on the boundary of the fundamental-mode dominance. Arrows indicate that 
fundamental- mode stability requires larger magnetic fields than those we simulated. In the region 
below the lines marked in long dashes (drawn only for illustrative purposes) the fundamental 
mode dominates. 
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This picture is qualitatively different from that drawn by Toth & Draine. For a very rough 
approximation to the interstellar-medium cooling function, they consider the simple broken power 
law shown in Figure |l|. Their linear analysis then yields for the fundamental-mode stability limit 
the line drawn in short dashes in Figure H- The results for our cooling function indicate that much 
larger magnetic fields are required to change the stability of the oscillations. The major difference 
between the cooling functions is that ours is much steeper (and hence more unstable) for shocks 
with Vs < 200kms~^, but flatter (hence more stable) for higher- velocity shocks. Our picture 
agrees with the results from single power-law cooling above: the fundamental mode dominates the 
oscillations, and the oscillations have significant amplitudes in the presence of typical magnetic 
fields when a < 0. 



3. Conclusions 

In our simulations the behavior of the overstability is determined by the typical shock 
temperature (hence the mean shock speed) and the upstream magnetic field. We divide the 
unstable cases into two categories: either the fundamental mode is unstable, or it is absent 
and the instability is dominated by the first overtone mode. The periods of these observed 
modes are roughly consistent with the prediction of linear perturbative analyses ( phevalier fc 



Imamura 1982| ; Toth fc Draine 19931 ). Simulations with cooling functions (eq. ||3|) that are 



truncated one-component power laws show a strong dependence on power-law index a and the 
Alfven Mach number M^^. The earlier simulations presented in Toth & Draine ( |1993| ) all 
have a > 0. For typical magnetic field strengths they found oscillations of limited (^ 10%) 
amplitudes dominated by overtone modes. For similar magnetic fields but a < 0, we observe that 
the fundamental mode dominates and that the oscillations have large amplitudes. We note a 
significant drop in oscillation amplitude accompanying the switch from fundamental to overtone 
mode as the magnetic field is raised, as we show for the typical case of a = —0.8. Even when 
the overtone mode dominates, the oscillations have substantial, nonlinear amplitudes. The linear 
analysis of Toth & Draine yields the values of the reciprocal Mach numbers M^^p and M^^q at 
which the fundamental and overtone modes become linearly stable for the full range of a. Our 
nonlinear simulations do not disagree with the linear analysis: the reciprocal Mach number 
at which the dominant oscillations change from fundamental to overtone mode differs significantly 
from M^^p, but always satisfies < M^^p. 

These single power-law cooling results explain our simulations that use a more realistic, 



piecewise power-law cooling function based on the results of Raymond et al. ( 1976 ). These 
simulations encompass mean shock speeds between 150 and 240kms~^ and Alfven velocities up 
to 20kms^^. Over most of the parameter space we observe large-amplitude fundamental-mode 
oscillations. The instability is so strong that quite substantial magnetic fields {va ^ lOkms"^) 
are required to change the qualitative behavior, and then only at the lowest and highest shock 
speeds. Again, we contrast our qualitative picture with the discussion of Toth & Draine, who 
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consider a very rough model for the interstehar-medium coohng function with a single power 
law (a = —1/2) above T = 10^'^ K. With a linear analysis they find that magnetic fields of 
realistic amplitudes may stabilize radiative shocks with velocities up to approximately 175 kms~^ 
in the warm ionized medium and the warm neutral medium. For some cases, we find that fields 
larger than any we simulated — that is, va > 40kms~^ — would be required to affect the stability 
behavior. The discrepancy between the two treatments points out that an accurate determination 
of the magnetic field required to stabilize shocks requires a more careful treatment of the cooling 



than the one we have attempted (see, for example, Innes 1992 ). However, the minimum magnetic 
field to stabilize the shock is quite dependent on the local power-law index a of the cooling 
function, and in any realistic treatment we expect a to vary over the velocity range of interest 
(130 kms~^ ^ (vs) <^ 250 kms~^). We then anticipate that the stabilizing magnetic field will vary 
from case to case but could be quite large, as found in our simulations. 

It must be noted that some aspects of gas behavior are not well modeled in our simulations. 
Perhaps the most prominent missing element is an accurate treatment of cooling, incorporating 
ionization and recombination in the cooling region, as well as in the upstream preshock gas and 
the cooled postshock gas. In particular, our cooling function only crudely mimics the turnoff of the 
cooling law upon recombination. There may as well be geometrical effects that require extending 
this work into two or three spatial dimensions. Further work will need to address one or more of 
these deficiencies. 
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Fig. 1. — The simulations use the coohng function indicated by the sohd hne, a piecewise power-law 
fit to the result of Raymond et al. ( 1976| ), and arranged to turn off a low temperatures. The dashed 
line shows the simple function adopted by Toth & Draine ( 1993 ) for application to the interstellar 
medium. The four dots indicate the postshock gas temperatures for shocks moving at 150, 180, 
210, and 240 km s"^ 
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Fig. 2. — From our simulations with power-law cooling, we find that the fundamental mode 
dominates the instability roughly in the region below the long dashed lines in (a, M^^) parameter 
space, where a is the power-law index. (The short line segments indicate simulations between 
which the fundamental mode becomes stable.) The short dotted lines indicate curves of marginal 
stability for the fundamental (Mj[^, lower) and overtone (M^q, upper) modes as calculated by 
Toth & Draine ( 1993 ). The dots indicate the twelve simulations discussed by Toth &; Draine. 
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Fig. 3. — Typical sizes of shock fluctuations (relative to steady-state values) as a function of 
reciprocal Alfven Mach number, for power-law cooling (eq. |^]) with a = —0.8. The fundamental 
mode dominates for M'^ < 0.08; and the overtone, for M^^ > 0.1. 
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Fig. 4. — Shock velocity (in units of lOOkms"^) as a function of time for mean shock speeds 
150, 180, 210, and 240kms~^ (from left to right) and upstream Alfven velocities 0, 2, 5, 10, and 
20kms~^ (from bottom to top). 
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Fig. 5. — From our simulations we find the fundamental mode is unstable roughly in the region 
below the long dashed lines. (The short line segments indicate simulations between which the 
fundamental mode becomes stable; the transition happens at magnetic fields above the points 
marked with arrows.) For comparison, the region below the short dashed lines is found to be 
unstable to the fundamental mode by Toth &; Draine ( 1993| ). The difference is attributable to the 
different cooling function (cf. fig. |l]). 
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Table 1: Average shock speeds in simulations 


Intended average 


Upstream Alfven velocity 


Inflow velocity Actual average 1 


(relative) 


( kms~^) 


(kms-^) 


(kms-^) 


( kms""^) 


150 





149 


150 ± 1 


150 


2 


148 


149 ± 1 


150 


5 


146 


151 ±3 


150 


10 


142 


150 ± 1 


150 


20 


134 


150 


180 





179 


181 ± 1 


180 


2 


178 


182 ±2 


180 


5 


176 


181 ± 1 


180 


10 


172 


180 ±2 


180 


20 


157 


177 ±2 


210 





209 


210 ± 1 


210 


2 


208 


210 ± 1 


210 


5 


206 


210 ± 1 


210 


10 


202 


211 ±3 


210 


20 


195 


210 ±2 


240 





239 


240 ± 1 


240 


2 


238 


238 ±2 


240 


5 


236 


238 ±3 


240 


10 


232 


240 ±1 


240 


20 


225 


240 ±2 
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Table 2: Typical parameters of shock-velocity oscillations without magnetic fields 



Mean shock velocity Velocity fiuctuation Fluctuation period^ Linear period Tp 
(kms^^) (kms"''") (yr) (yr) 



130 


80 


4 400 


4 000 


140 


95 


7 000 


4 500 


150 


140 


11 000 


5 600 


160 


140 


13 000 


7 900 


170 


160 


17 000 


11 000 


180 


160 


20 000 


16 000 


190 


160 


30 000 


22 000 


200 


160 


34 000 


29 000 


210 


200 


40 000 


35 000 


225 


210 


47 000 


43 000 


240 


240 


63 000 


53 000 



"nif = 1 cm ^; all times vary as ri 



